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I. INTRODUCTION 



Many complex processes that contain multi-jet final states will be observed at the Large 
Hadron Collider (LHC). Detailed studies of these processes will give us information about 
the mechanism of electroweak symmetry breaking and, hopefully, reveal physics beyond the 
Standard Model. To assess this information and interpret experimental data correctly, an 
accurate theoretical description of processes at the LHC is required. In principle, Monte- 
Carloprograms based on leading order (LO) computations do provide such a description [l|, 
y, y, I4, 5]. However, as experience with the TEVATRON, LEP and HERA data has shown, 
LO predictions often give only rough estimates. To extract maximal information from data, 
more precise predictions based on next-to-leading order (NLO) calculations are required. 
NLO predictions are also instrumental for reliable estimates of theoretical uncertainties 
related to the truncation of the perturbative expansion. 

For sufficiently complicated final states, computations required for the LHC are very chal- 
lenging. In the standard approach one uses perturbative expansion of scattering amplitudes 
in terms of Feynman diagrams. This gives an algorithm suitable for numerical implementa- 
tion. However, even for tree amplitudes, the number of Feynman diagrams grows faster than 
factorial with the number of external particles involved in a scattering process. As a conse- 
quence, computing time needed to evaluate a scattering amplitude at a single phase-space 
point, grows at least as fast. This means that computational algorithms based on expansion 
in Feynman diagrams are necessarily of exponential complexity, an undesirable feature. 

In tree-level calculations exponential growth in complexity is avoided by employing re- 
cursion relations. These relations re-use recurring groups of off-shell Feynman graphs in an 
optimal manner js, y, y, y, 1^. The use of recursion relations for tree amplitudes leads 
to a computational algorithm of polynomial complexity so that computing time grows as 
some power of the number of external legs. Because of that, the problem of evaluating tree 
amplitudes is considered a solved problem by the high-energy physics community. As was 
pointed out in Ref. a polynomial-complexity algorithm is not available for one- loop 
computations; constructing such an algorithm is the goal of the present paper. 

When performing NLO computations in the Standard Model, many difficulties arise. We 
need to calculate both virtual one-loop corrections and real emission processes with one ad- 
ditional particle in the final state. Currently, the bottleneck in NLO calculations for multi- 
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particle processes is the computation of virtual corrections. The difficulty related to the 
factorial growth in the number of Feynman diagrams is further amplified by a large number 
of terms generated when tensor loop integrals are reduced to scalar integrals. Nevertheless, 
thanks to modern computational resources, standard methods based on Feynman diagrams 
and tensor integrals reduction [l^ may be extended brute force to deal with multi-particle 
processes. Striking examples of the success of this approach are recent computations of elec- 
troweak corrections to e~^e~ four fermions process [l^ and one-loop six-gluon scattering 
amplitudes Ref. Note, however, that a single phase space point for six-gluon scattering 
is evaluated in about nine seconds, which is 10,000 times slower than the evaluation time for 
four-gluon scattering amplitudes generated using the same procedure. It is clear that fur- 
ther application of brute force approaches to yet higher multiplicity processes are becoming 
unfeasible. 

Unitarity-based methods for multi-loop calculations were suggested as an alternative to 
the expansion in Feynman diagrams long ago. We review the status of these calculations 
in the next section. The goal of the present paper is to describe a polynomial-complexity 
computational algorithm for one-loop amplitudes that provides both cut-constructible and 
rational parts. We start with the idea of generalized unitarity in D-dimensions and de- 
velop it to an algorithm amenable to numerical implementation. Since one-loop amplitudes 
are built up from tree amplitudes, albeit in higher- dimensional space-time, the polynomial 
complexity of the algorithm is ensured. The method is flexible and can be applied to scat- 
tering amplitudes with arbitrary internal and external particles. In particular, dealing with 
massive particles is straightforward. 

The outline of the paper is as follows. In section II we give an overview of the current sta- 
tus of unitarity techniques. The structure of one-loop scattering amplitudes in D-dimensions 
is discussed in section III. Section IV contains the discussion of the D-dimensional residues 
and algebraic extraction of the coefficients of master integrals. In section V we apply the 
formalism to gluon scattering amplitudes. Numerical results for four-, five- and six-gluon 
scattering amphtudes are reported in section VI. We conclude in section VII. 
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II. THE STATUS OF UNITARITY METHODS 



Unitarity-based methods for loop calculations were suggested as an alternative to the 
Feynman-diagrammatic expansion long ago 1^, Q, [l^] • It was argued that for gauge theories 



these methods lead to higher computational efficiency than traditional methods [18|, Il9 |. 
Within unitarity-based methods, computations employ tree scattering amplitudes, rather 
than Feynman diagrams, thereby avoiding many complications. Unitarity cuts factorize one- 
loop amplitudes into products of tree amplitudes. Therefore, in numerical implementations 
computing time grows with the number of unitarity cuts, rather than Feynman diagrams, 
and depends on the efficiency of algorithms employed for evaluating tree amplitudes. 

The unitarity-based approach gives a complete description of one-loop scattering ampli- 
tudes if it is applied in D dimensions. Any dimensionally regulated multi-loop Feynman 
integral is fully reconstructible from unitarity cuts 17j. Clearly, in this case we have to 
associate D-dimensional momenta and polarization vectors with each cut on-shell line, to 
obtain one- loop amplitudes [2^. 

The first successful application of unitarity-based techniques for one-loop computations 
employed a four-dimensional variant of the unitarity-based approach, where four- dimensional 
states were associated with each cut line 2]|. In analytic calculations, such a procedure has 



the advantage of allowing full use of the spinor-helicity formalism. In this way, however, 
we obtain only the so-called cut-constructible part of the full one-loop amplitude. The 
missing part is referred to as the rational part and has to be determined by other methods. 
In particular, in supersymmetric theories, rational parts are known to vanish. In other 
cases, rational parts can be fixed by using factorization properties of one-loop amplitudes in 
coUinear limits 2l| . 



An important step in developing unitarity-based methods was made in Ref. 22| where the 
idea of generalized unitarity was introduced. In particular, it was shown that coefficients 
of four-point scalar integrals for multi-gluon processes can be calculated using quadruple 
cuts. This approach is also suitable for direct numerical implementation. A quadruple cut 
factorizes an one-loop amplitude into a product of four tree amplitudes. From unitarity 
constraints, we derive two complex solutions for the loop-momentum. The product of tree 
amplitudes can be evaluated using these solutions. Coefficients of scalar four-point functions 
are obtained by taking averages. 



4 



Unfortunately, the simplicity of the above procedure does not generalize easily to the 
computation of full one-loop amplitudes. In that case, we also have to determine coefficients 
of one- , two- and three-point scalar integrals. For example, when a double cut is applied to 
determine coefficients of two-point functions, we have to account for the fact that parts of 
these contributions are already contained in quadruple and triple cuts. Analytic separation 
of these overlapping contributions proved to be complicated. 

An efficient algebraic method to separate the overlapping contributions was outlined in 
Ref. [23|. Developing upon this idea, flexible unitarity-based computational techniques were 



suggested in Refs. 2J, |25[ where it was shown how generalized unitarit y ca n be implemented 
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27|. All these papers 



numerically. Analytic techniques were further developed in Refs. 
address primarily cut-constructible parts of scattering amplitudes. 

These developments provide a polynomial-complexity computational algorithm for cut- 
constructible parts of scattering amplitudes. However, techniques for calculating rational 
parts are much less developed. Three methods for computing rational parts are currently 
available. 

In Refs. 



291] it was suggested to determine rational parts of tensor one-loop integrals 
analytically and employ traditional diagrammatic methods to obtain scattering amplitudes. 
This leads to an exponential-, rather than polynomial-complexity algorithm, negating all 
the progress achieved with the determination of the cut-constructible part using numerical 
unitarity techniques. 

The other two methods are so far analytic, but both should in principle be suitable for 
numerical implementation. The so-called bootstrap method sets up a recursive computation 
of rational parts \3l\ similar to tree- level unitarity-based recursion relations [10|. How- 
ever, in its current formulation the bootstrap approach is not directly suitable for numerical 
implementation since both the cut-constructible and rational parts contain spurious poles. 
When the two parts are added together, spurious poles cancel. Unitarity-based recursion 
relations for the rational part can only be constructed once spurious poles are removed from 
the rational part. The procedure to remove these poles from the rational part is called 
cut-completion; it requires analytic knowledge of the cut-constructible part. Clearly, for 
complicated multi-particle processes this is a serious disadvantage. 



The other approach is a variant of D-dimensional unitarity method 
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33| . Since the cut 



lines are D- dimensional, tree scattering amplitudes have to be calculated in D-dimensions. 
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FIG. 1: The generic A^-point loop amplitude. 

Currently, this method is restricted to analytic applications for purely gluonic one-loop scat- 
tering amplitudes with a scalar particle propagating in the loop. In the following sections we 
describe a computational method based on /^-dimensional unitarity that permits straightfor- 
ward calculation of both cut-constructible and rational parts of arbitrary one-loop scattering 
amplitudes. 



III. ONE-LOOP AMPLITUDES AND DIMENSIONALITY OF SPACE-TIME 



Since one-loop calculations in quantum field theory lead to divergent expressions, we 
require regularization at intermediate stages of the calculation. Such regularization is ac- 
complished by continuing momenta and polarization vectors of unobserved virtual particles 
to D ^ A dimensions 3J|. On the other hand, it is convenient to keep momenta and po- 
larization vectors of all external particles in four dimensions since this allows us to define 
one-loop scattering amplitudes through helicities of external particles. Once the dependence 
of a one-loop amplitude on the dimensionality of space-time is established, we interpolate 
to a non-integer number of dimensions D = 4 — 2e. The divergences of one-loop amplitudes 
are regularized by the parameter e. 

We wish to arrive at a numerical implementation of this procedure. To this end, it is 
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crucial to keep the number of dimensions in which virtual unobserved particles are allowed 
to propagate as integer since only in this case loop momenta and polarization sums of 
unobserved particles are fully defined. Therefore, we determine the dependence of one-loop 
amplitudes on the dimensionality of space-time treating the latter as integer and arrive at 
non-integer values (e.g. D = A — 2e) later on, by simple polynomial interpolation. 

Any D-dimensional cyclic-ordered A^-particle one-loop scattering amplitude (Fig. [1]) can 
be written as 

^^^^^ - J zRW^^^d^^TTd^ ' 
where {pi} and { Jj} are the two sets that represent momenta and sources (polarization vec- 
tors, spinors, etc.) of external particles. The numerator structure J^{{pi}, {Ji}',l) depends 
on the particle content of the theory. The denominator is a product of inverse propagators 

di = di{l) = {l + qif -m'^i = ^-qo + Y^Pi^ - , (2) 

where the four-vector go represents the arbitrary parameterization choice of the loop mo- 
mentum. 

The one-loop amplitude can be written as a linear combination of master integrals 

^ ^^/^ '-'n«2«3-'jii2j3 ^ / J "n»2-'ni2 ^ / ^ "n-'n ; \^) 

where we introduced the short-hand notation = 1 < ii < 12 < ■ ■ ■ < in N . The 

master integrals on the r.h.s. of Eq. ([3]) are defined as 

/(^). = I J-1 i (4) 

J z(7r)^/2c/,,---c/,,, • ^ ^ 

The coefficients of master integrals for this choice of basis depend on the number of dimen- 
sions D which, in practical calculations in dimensional regularization, needs to be taken as 
D = A — 2e. Since we aim at numerical implementation of /^-dimensional unitarity, this 
is inconvenient. We explain below how to change the basis of master integrals to make 
coefficients D-independent. 

The D-dependence of one-loop scattering amplitudes associated with virtual particles 
comes from two sources. When we continue loop momenta and polarization vectors to higher- 



dimensional space time, the number of spin eigenstates changes. For example, massless spin- 
one particles in Dg dimensions have Dg — 2 spin eigenstates while spinors in Dg dimensions 
have 2(^=-2)/2 

spin eigenstates. In the latter case, Dg should be even. 
The spin density matrix for a massless spin-one particle with momentum I and polariza- 
tion vectors efi' is given by 



E^?(0e«(0 = -^7i^=^ + ^^^^^ (5) 

i=l 

where is an arbitrary light-cone gauge vector associated with a particular choice of po- 
larization vectors. Similarly, the spin density matrix for a fermion with momentum / and 
mass m is given by 

2(Ds-2)/2 ^ 

E u« = /+m = E/^7^ + m . (6) 

While, as we see from these examples, the number of spin eigenstates depends explicitly on 
the space-time dimensionality, the loop-momentum / itself has implicit D-dependence. We 



can define the loop momentum as a D-dimensional vector, with the requirement D < Dg |35|. 
We now extend the notion of dimensional dependence of the one-loop scattering amplitude 
in Eq. ([T]) by taking the sources of all unobserved particles in Ds-dimensional space-time 

AinMiP^}^ W) - J d,d,...d^ ■ 

The numerator function N'^^''\{pi]^{Ji];l) depends explicitly on Dg through the number 
of spin eigenstates of virtual particles. However, the dependence of the numerator function 
on the loop momentum dimensionality D emerges in a peculiar way. Since external particles 
are kept in four dimensions, the dependence of the numerator function on — 4 components 
of the loop momentum / appears only through its dependence on P . Specifically 

i'=f-T' = il-il-il-il-Y^ih (8) 

where I and / denote four- and (D — 4)-dimensional components of the vector /. It is apparent 
from Eq. ([8]) that there is no preferred direction in the {D — 4)-dimensional subspace of the 
D-dimensional loop momentum space. 

A simple, but important observation is that in one-loop calculations, the dependence of 
scattering amplitudes on Dg is linear. This happens because, for such dependence to appear, 
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we need to have a closed loop of contracted metric tensors and/or Dirac matrices coming 
from vertices and propagators. Since only a single loop can appear in one-loop calculations, 
we find 

Af(^^\l)=Af,{l) + {D,-4)m). (9) 

We emphasize that there is no explicit dependence on either Ds or D in functions A/o,i. 

For numerical calculations we need to separate the two functions A/q,!. To do so, we 
compute the left hand side of Eq. ([9]) for Dg = Di and Dg = D2 and, after taking appropriate 
linear combinations, obtain 

(A - 4)Ar(^i)(/) - [Di - 4)Ar(^2)(/) 



D2-D1 



Because both Di and D2 are integers, amplitudes are numerically well-defined. We will 
comment more on possible choices of D12 in the forthcoming sections; here suffice it to say 
that if fermions are present in the loop, we have to choose even Di and D2. 

Having established the Dg-dependence of the amplitude, we discuss analytic continuation 
for sources of unobserved particles. We can interpolate Dg either to Dg 4— 2e (the t'Hooft- 
Veltman (HV) scheme) 3J] or to — 4 (the four-dimensional helicity (FDH) scheme) 35|. 
The latter scheme is of particular interest in supersymmetric (SUSY) calculations since all 
SUSY Ward identities are preserved. We see from Eq. ([9]) that the difference between the 
two schemes is simply — 2eA/i. 

We now substitute Eq.([Tni) into Eq. (JTj). Upon doing so, we obtain explicit expressions 
for one-loop amplitudes in HV and FDH schemes. We derive 

= - {d^) {Ad,d.^d.) - A,o,D.^o.)) . (11) 

We emphasize that Dg = D12 amplitudes on the r.h.s. of Eq. (fTTi) are conventional one-loop 
scattering amplitudes whose numerator functions are computed in higher-dimensional space- 
time, i.e. all internal metric tensors and Dirac gamma matrices are in integer Dg = Di^2 
dimensions. The loop integration is in D < Dg dimensions. It is important that explicit de- 
pendence on the regularization parameter e = {4 — D)/2 is not present in these amplitudes. 
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For this reason, Eq. (fTTi) renders itself to straightforward numerical implementation. In par- 
ticular, numerical implementation of Ds-dimensional unitarity cuts is now straightforward, 
as cut internal lines possess well-defined spin density matrices. The D^-dimensional unitar- 
ity cuts, that we discuss in detail in the next section, decompose amplitudes into a linear 
combination of master integrals, see Eq. We can choose the basis of master integrals 
in such a way that no explicit D-dependence in the coefficients appears. Only after the re- 
duction to master integrals is established, we continue the space-time dimension associated 
with the loop momentum to — 4 — 2e, thereby regularizing the master integrals. 



IV. ^^-DIMENSIONAL UNITARITY CUTS 



The amplitudes on the r.h.s. of Eq. (ITTI) are most efficiently calculated by using gen- 
eralized unitarity applied in Dg dimensions. Since many aspects of the calculation in this 
case are the same as in four-dimensional generalized unitarity, we focus on new features that 



appear for Ds > 4. Our discussion follows Ref. [2J] which details the = 4 case. 

The important issue is the set of master integrals that we have to deal with when applying 
i^s-dimensional unitarity cuts. Cutting a propagator i requires finding the loop momentum 
/ such that equation di{l) = is satisfied. Therefore, each cut imposes one constraint on the 
loop momentum. In four dimensions, where loop momentum has only four components, one 
can cut at most four propagators without over-constraining the system of equations. This 
leads to the conclusion that four-point integrand functions in four-dimensions are required 
whereas five-point and higher-point functions aren't. 

When we compute an amplitude in > 4 dimensions, we may also use the additional D— 
4 components of the loop momenta to set more inverse propagators to zero. However, since 
all external momenta are four-dimensional, additional components of the loop momentum 
enter all propagators in a particular combination 



D D 

2 



1=5 i=5 

where nj>4 are orthonormal basis vectors {ni-nj = 6ij) of the — 4)-dimensional sub-space 
embedded in a D-dimensional space. Therefore, when we move from four to Dg dimensions, 
at most five inverse propagators can be set to zero. Therefore, for Dg > 4, the highest- 
point master integral that should be included into the master integrals basis is a five-point 
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function. Hence the integrand of the A^-particle amphtude in Eq. ([T]) can be parameterized 



as 



2-^ A- A. A- rl- ^. /7. W 



did2 ■ ■ ■ d]\[ di^di^di^di^di^ , di-^di^di^di^ 

where the dependence on the external momenta and sources are suppressed. From four- 
dimensional unitarity we know that computation of each cut of the scattering amplitude is 
simplified if convenient parameterization of the residue is chosen. We now discuss how these 
parameterizations change when Dg-dimensional unitarity cuts are considered. 



A. Pentuple residue 

To calculate the pentuple residue, we choose momentum / such that five inverse propa- 
gators in Eq. f|T3l) vanish. We define 

The momentum kjkmn satisfies the following set of equations di{lijkmn) = ■ ■ ■ = dnilijkmn) = 
0. The solution is given by 

where m„ is the mass in the propagator dn which is chosen to be as dn = by adjusting 

the reference vector go- The parameters an can be chosen freely. The four-dimensional vector 
depends only on external momenta and propagator masses. It is explicitly constructed 



using the Vermaseren-van Neerven basis as outlined in Ref. [2J]. The D — A components 
of the vector UjUmn are necessarily non-vanishing; for simplicity we may choose lijkmn to be 
five-dimensional, independent of Dg. We will see below that this is sufficient to determine 
pentuple residue. 

To restrict the functional form of the pentuple residue e-ijUmni}) we apply the same reason- 
ing as in four- dimensional unitarity case, supplemented with the requirement that e^^kmni^) 
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depends only on even powers of Se] this requirement is a necessary consequence of the discus- 
sion around Eq. ([8]). These considerations lead to the conclusion that the pentuple residue 
is independent of the loop momentum 

(16) 

To calculate e*^°^ in the FDH scheme, we employ Eq. (1101) and obtain 

(o),FDH _ ^ ^2 - 4 ^ ^ f^{^!^\_fJ^i:i±\T, f^n^\ 

The calculation of the residues of the amplitude on the r.h.s. of Eq. f[T7|) . is simplified by 
their factorization into products of tree amplitudes 

ReSijkmn ( ^ ^ j = 2j-^(^i;P*+l, • • -yPj, "^j) X M{lj;Pj+l, . . . , Pk] -Ik) (18) 

xM{lk]Pk+l, ■ ■ ■,Prn] -Im) X M{lm\Pm+l, ■ ■ ■ , Pn, -Q X M^ln] Pn+1- ■ , Ph -k) ■ 

Here, the summation is over all different quantum numbers of the cut lines. In particular, we 
have to sum over polarization vectors of the cut lines. This generates explicit Dg dependence 
of the residue, as described in the previous section. Note that the complex momenta = 
+ (Ih are on-shell due to the unitarity constraint dh = 0. 



B. Quadrupole residue 

The construction of the quadrupole residue follows the discussion of the previous subsec- 



tion and generalizes the four- dimensional case studied in 23|, |2J] . We define 



Ar(^^)(/) 



C (0 - Res,,,. _ _ _ d^A^AA. j ' ^^^^ 

where the last term in the r.h.s. is the necessary subtraction of the pentuple cut contribution. 
We now specialize to the FDH scheme. In this case, the most general parameterization of 
the quadrupole cut is given by 

where si = / ■ ni. We used the fact that, in renormalizable quantum field theories, the 
highest rank of a tensor integral that may contribute to a quadrupole residue is four and 
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that only even powers of Se can appear on the r.h.s of Eq. (120|) . The solution of the unitarity 
constraint is given by 



D 



h=5 



The vector V4 is defined in the space spanned by the three independent inflow momenta 

jhe one-dimensional "transverse space" , 



24|. 



{ki, kj, kk} and rii is the unit vector that describes 
i.e. ki ■ ni = 0, kj ■ ni = 0, kk ■ ni = 0, rii ■ rij = Sij 

To determine the momentum- independent coefficients df^'i!'^"'^^ in Eq. (j20l) . we compute 
the quadrupole residue of the one-loop scattering amplitude for different values of the loop 
momentum / that satisfies the unitarity constraint. This entails choosing different values for 
parameters in Eq. (12T1) . As a first step, we may choose / to be a four-dimensional vector 
embedded in D^-dimensional space (aj>5 = 0). Then Se = and the parameterization of the 
becomes identical to a four-dimensional case. Standard manipulations 



residue in Eq. ((20 
described in Refs. 



23 



24( 1 then allow us to find d^^^l^ and c?-J^^. To determine the remaining 



coefficients, we consider loop momenta in dimensions D > 4 such that di = dj = dk = dm = 
but Se 7^ 0. We accomplish this by adjusting the value of a^. By choosing appropriate loop 
momenta, we determine all the remaining three coefficients of the quadrupole residue by 
solving a linear system of equations. 



C. Triple , double and single-line cuts 

Calculation of triple, double and single cuts proceeds in full analogy with what has been 
described for pentuple and quadrupole cuts. The only modification concerns parameteriza- 
tion of residues. 

The general parameterization of a triple cut in the FDH scheme is given by 

-FDH/n _ (0) (1) (2) (3) . 2 „2n , / (4) (5) (6) x 

Cijk ['■) - '^ijk + '^ijk^i + '^ijk^2 + C^ji^[S;i^ - S2) + SiS2[C^j,^ + C^ji^Si + C^j,^S2) 

I „(7) „ „2 , (8) 2 , (9) 2 /99\ 

where Si = I ■ ni and S2 = I ■ The solution of the unitarity constraint is given by 

'5. = v-," + („,„,. + + f . (23) 



h=5 
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The vector V3 is defined in the space spanned by the two independent inflow momenta {ki, kj} 
and ni^2 are orthonormal vectors that describe the two-dimensional "transverse space", i.e. 



ki ■ nio = 0, kj ■ ni 9 = 0, ■ n,- 



■'V 



The general parameterization of a double cut is given by 



4) + bf'^S2S, 



(24) 



where Si = I ■ ni, S2 = I ■ n2 and S3 = / ■ 71.3. The solution of the unitarity constraint is given 
by 

D 



I' = vr + 



^2 



m2 



0^2 + ^3 + tts 



(25) 



/i=5 



The vector V2 is proportional to the inflow momentum {ki} and 71-1.2,3 are orthonormal vectors 



that describe the three-dimensional "transverse space", i.e. k^ ■ 'n.1,2,3 = 0, ■ rij = dij 



24|. 



Finally, we note that the parameterization of a single-line residue is the same as in the 
four-dimensional case described in Ref. 24 1. 



D. One-loop amplitudes and master integrals 

Following the strategy outlined in previous subsections, we obtain pentuple, quadrupole, 
triple, double and single-line residues. These residues give us coefficients of master integrals 
through which the one-loop amplitude is expressed. 

Because we have more coefficients in our parameterization of residues, than in the four- 
dimensional case, we end up with a larger number of master integrals. The new master 
integrals include the five-point function and various Sg-dependent terms that appear in 
quadrupole, triple and double residues. We will comment on the fate of the five-point 
function shortly. 

Consider now the Sg-dependent master integrals. Some of them contain scalar products 
/ ■ ni^2 and vanish upon angular integration over /. Neglecting these spurious terms, we 
have to deal with four additional master integrals. We can rewrite those integrals through 
conventional four-, three- and two-point functions in higher-dimensional space-time. We 
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find 



f d^l sj ^ {D^2){D-A) 

J {t7i)D/2 d,^d,A^d^, 4 ^ ^ 

J {t7r)D/^ d,AAs 2 
/• d^/ _ (/^-4) ^^, 

7 {tTr)D/2diA2 2 • 

Using Eq. (l26l) . we arrive at the following representation of the scattering amplitude 

A. = V p(0) jiD) 

[«l|«5] 

+ V r/(°^ y^^) - ^ ~^ r/^^^ y^^^^^ + {D-4:){D - 2) (4) (^^+4) 

' / J 1 ili2i3i4 411223*4 2 41*24344 41424344 /| "41424344 -'41424344 

[4l|44] ^ 

+ V l^r(°) y^^) - ^ ~ ^ r^^) /-(^+2)^ 

/ ^ I Wl4243 414243 414243 414243 I 



[41I44] 

[ilK3] ^ 

+ E 4S - ASr') + E «f 4f' . (27) 

[4l|42] ^ ^ 41 = 1 

We emphasize that the explicit D-dependence on the r.h.s. of Eq. ( !27|) is the consequence 
of our choice of the basis for master integrals in Eq. fl26|) . 

We note that the above decomposition is valid for any value of D. We can now inter- 
polate the loop integration dimension D to — 4 — 2e. The extended basis of master 
integrals that we employ provides a clear separation between cut-constructible and rational 
parts of the amplitude. The cut-constructible part is given by the integrals in D-dimensions 
in Eq. fl27j) . while the rational part is given by the integrals in D + 2 and D + 4 dimen- 
sions. However, it is possible to use smaller basis of master integrals by rewriting integrals 

{4S^X^lS24^^S^3^^^hi2^^} in *erms of {I-SLiA^ ^SL^ ^Sl} ™g integration-by- 
parts techniques. 

Since we are interested in NLO computations, we only need to consider the limit e 
in Eq. fl27|) and neglect contributions of order e. This leads to certain simplifications. First, 
in this limit, we can re- write the scalar 5-point master integral as a linear combination of 
four-point master integrals up to 0{e) terms. If we employ this fact in Eq. fl271) . we obtain 

e™..,, 4°'. + E 4°'..) = E + 0(«) . (28) 

[41I44] y [41I44] 
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Note that since scalar five-point function is cut-constructible, it does not contribute to 
tlie rational part. Second, since lim£)-+4(-D — 4) x I^^lj^ = 0, tlie 6-dimensional 4-point 
integral can be neglected. Finally, the 6-dimensional 2-point, 6-dimensional 3-point and 8- 
dimensional 4-point integrals are all ultraviolet divergent and produce contributions of order 
e~^. Therefore, factors {D — 4) that multiply these integrals in Eq. flTTI) pick up divergent 
terms and produce e-independent, finite contributions 

{D-4){D-2) (^+4) ^ _1 

4 ^n^2^3^4, 3' 



(£-4) , . 

lim = ^ + - (g, - g,) . 

Combining everything together, we can write any one-loop amplitude up to terms of order 
e as a linear combination of cut-constructible and rational parts. 

An = A%'^ + Rn. (30) 
The expression for the cut-constructible part reads 

N 

jCC ^ V- T(0) Ai-2e) V- (0) .(4-2.) ,(0) .(4-2.) , „(0) j(4-2e) .^iN 

[«i|«4] [ii\is] [hlh] il=l 

where coefficients df^l.i^ are implicitly defined in Eq. fl28p . For the rational part, we obtain 

3 2 Y 6 2 j 

[n|«4] [iilia] [n|«2] 

Finally, we point out that various master integrals required for one-loop computations can 



be found in Ref. [36 1. 



V. GLUON SCATTERING AMPLITUDES IN QCD 

We applied the method described in the previous sections to compute gluon scattering 
amplitudes in QCD. It is well-known that these amplitudes can be represented as linear 
combinations of simpler objects, called color-ordered sub-amplitudes. For example, the tree 



amplitude for n-gluon scattering reads 



37 



38 
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where pi,Xi,ai stand for momenta, helicities and color indices of external gluons, g is 
the coupling constant and are generators of the SU(Ai'c) color algebra normalized as 
Tt{T"'T^) = 6°'''. The sum in Eq. fl33l) runs over (n — 1)! non-cyclic permutations of the set 
{1, . . . ,n}. Amplitudes A^^'^'^ {ptli)' ^ ■ ■ ■ ^ptln)^ color-ordered tree sub-amplitudes. 

One-loop amplitudes can be decomposed in a similar, but slightly more complicated 
fashion 39|]. Considering one- loop amplitudes for n-gluon scattering, where all internal 



particles are also gluons, we write 



where 



and 



[n/2]+l 

An{{Pi,Xi,ai}) = g'' Yl Yl Gr„.c((T)A,c(o-), (34) 



Gr^.i (a) = NJTi (T^^w ■ ■ ■ T'^-t") ) , (35) 



Gr„.e(a) = Tr {T^-w ■ ■ ■ T'^-c-d ) Tr {T^-(^^ ■ ■ ■ T''-(") ) . (36) 

In Eq. f lM|) [x] is the largest integer number smaller than or equal to x and Sn.c are subsets 
of Sn, that leave the double trace in Eq. (!36l) invariant. 

Since sub- amplitudes possess a number of symmetries under cyclic permutations of ex- 
ternal particles and parity, not all An.i amplitudes are independent. Moreover, all An,c>i 
amplitudes can be written as linear combinations of An,i amplitudes. Because of that, there 
are four independent helicity sub-amplitudes for four- and five-gluon scattering and eight 
independent helicity sub-amplitudes for six-gluon scattering. Scattering amplitudes for four, 
five and six gluons are available in the literature; this allows us to check the validity of our 
method. 

The method described in the previous section is amenable to straightforward numerical 
implementation. We construct a list of all possible cuts of a given sub-amplitude. Those 
cuts that correspond to two-point functions for light-like incoming momenta are discarded 
since corresponding master integrals vanish in dimensional regularization. 

Each cut in the list is computed as a product of tree amplitudes which are obtained using 
Berends-Giele recurrence relations The recurrence relations themselves do not change 
when we construct gluon scattering amplitudes in higher-dimensional space-time. However, 
polarization vectors for cut gluon lines have to be extended. To discuss this extension 
explicitly, we now choose specific values for space-time dimensionalities D12 in Eq. (fTTl) . 
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Since we deal with pure gluonic amphtudes, we may consider Di = 5 and D2 = 6. We 
obtain 

^FDH ^ 2A[D,Ds=5) - -^{0,0^=6)- (37) 

Computation of residues discussed in the previous section requires Sg 7^ 0. This can be 
easily accomplished by allowing the projection of the loop momentum on the fifth direction 
to be non-vanishing, while always keeping I ■ uq = 0, even for Ds = 6. Hence, for computing 
residues of gluon amplitudes we have to consider a few cases of how four- and five- dimensional 
loop momentum can be embedded into five- and six- dimensional space-time. For the sake 
of clarity, we describe those cases separately. 

For Dg = 5 and four-dimensional loop momentum, we have I ■ = 0. This allows us to 
write 

p^^'i^.r]), /i, 1/ e 4 dim; 



E 

1=1 



(38) 



otherwise, 



where 



-9 



I ■ r] 



and ?7 is a four-dimensional light-cone vector such that / ■ 77 7^ 0. 

For Dg = 6 and four-dimensional loop momentum, we have / ■ 
choose the following expression for gluon density matrix 



(39) 



I ■ riQ = 0. We then 



i=l 



p^^{l, 77), /i, G 4 dim; 

-n'i^nl, /i = 5, = 5; 

-n(^<, = 6, 1/ = 6; 

0, otherwise. 



(40) 



For Z^s = 5 and five-dimensional loop momentum, we write = + Pn'^, where l-n^ = 0. 
Also, = 0, but / = — 7^ 0. Then, we construct polarization vectors by taking them to 
be in the three-dimensional subspace of a five-dimensional space, which is orthogonal to I 
and 77.5. Specifically, we have 



i=l 



-ujf"'{l); /i, z/ e 4 dim; 
0, otherwise, 



(41) 
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where to^'^ = —g^" + I'^f/f^ with all indices restricted to four-dimensions. From a four- 
dimensional viewpoint, Eq. (l4Ti) is a density matrix of a spin-one particle with the mass 
f. 



For Ds = Q and five-dimensional loop momentum, we take I = f + Pn'^- In this case 
/ ■ ng = 0. With this choice, all we need to do is to add one more polarization direction to 
Eq. dUD. We obtain 

f _ 

-tu^^'il); /i, G 4dim; 

-«, ;i = 6,z/ = 6, (42) 
0, otherwise. 

We now make the following observation. For a chosen loop momentum /, be it four- or 
five-dimensional, the calculation in Ds = Q differs from the calculation in = 5 by a single 
polarization component of the gluon, denoted ng in Eqs. (I40f42p . All external momenta 
and polarizations are four-dimensional, and the cut loop momentum / that satisfies the 
unitarity constraint is at most five- dimensional. Because of this the uq polarization gives 
non-vanishing contribution when it is contracted with itself through a metric tensor. From 
this point of view, its contribution is equivalent to that of an additional (real) scalar particle 
in the loop. To see this more explicitly, consider a three-gluon vertex with two gluons 
polarized along the sixth dimension in six- dimensional space and the third gluon with the 
four-dimensional polarization vector. Since none of the gluon momenta have a component 
along the sixth dimension, we obtain 

VSlJk^^ ^2, fc3)ernf ~ e^,\h - h^- (43) 

The object on the r.h.s. of the above equation is the scalar-scalar-gluon vertex. Similar 
considerations applied to a four-gluon vertex with two gluons polarized along the uq direction 
immediately lead us to conclude that it becomes a scalar-scalar-gluon-gluon vertex. Hence, 
any tree n-gluon amplitude with two gluons whose polarization vectors are = is 
equivalent to a tree scattering amplitude of n — 2 gluons and two scalars. This allows us to 
write 

where the amplitude Af) is computed with a particular propagator for all internal gluon 
lines —171^71^/1'^ which, as we argued above, is equivalent to scalar contribution. 



i=l 



U V 
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If we use Eq. ( H^ in Eq. (1371) . we derive the following result 

^™" = Aad.=5)-4i»- (45) 

It tells us that, from the point of view of generalized unitarity, computations of gluon 
scattering amplitudes in the four-dimensional helicity scheme are equivalent to calculations 
of gluon scattering amplitudes in five- dimensional space-time up to a term that can be 
associated with a scalar contribution. 

Finally, we may now write explicitly the relation between amplitudes computed in FDH 
and HV schemes. Using Eq. ( ITTi) and Eq. (jB]), we derive 

A'''' = ^(D,D.=5) - (1 + 26)4^) . (46) 

We see that the difference between the HV and FDH schemes is entirely due to the addi- 



tional "scalar" degree o 
re-derive the results of 



freedom that contributes to one-loop amplitudes. Note that we 



lol without resorting to SUSY arguments. The relation between 
unrenormalized amplitudes in two schemes becomes 

^FDH _ ^HV ^ 2e^^ (47) 
The contribution of a real scalar to one-loop gluon scattering amplitudes reads 

= ^Arec + finite , (48) 
be 

where 

r(i + .)r(i-.)V 

^ (4^)2- ■ ^'^J 

is the usual normalization factor and /i is the scale that maintains the dimensionality of loop 
integrals after the loop momentum is continued to -D = 4 — 2e dimensions. Neglecting terms 
of order e, we obtain the well-known relation between gluon amplitudes computed in FDH 



and HV schemes 40 1 



VI. RESULTS 



We now present the results of the calculation of four-, five- and six-gluon scattering 
amplitudes in QCD. If, for a given choice of gluon helicities, the corresponding tree amplitude 
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does not vanish, the n-gluon one-loop sub-amphtude is written in the following way 

=Cr (^-^ + 1 (^-y + + AA„...,A„j Arec(l^\...,n^"), 

(51) 

where Sj^j+i = 2pi ■ pi^i + i6, Pn+i = Pi- The convenience of representing scattering am- 
plitudes as in Eq. fl5T|) is that phase conventions for gluon polarization vectors drop out 
when functions Aai,,,,,a„ are computed; this feature allows for direct comparison with the 
literature. 

We have verified that our calculations correctly reproduce the divergent parts of Eq. (15T|) . 
In Tables [H |TT] and IIIII we give the results for finite parts of color-ordered sub-amplitudes 
Aai,...,a„- The finite parts are split into cut-constructible and rational parts 

Aai,...,A„ = ^T,...,x,, + 

for kinematic points specified below. 

On the other hand, some color-ordered sub-amplitudes for n-gluon scattering are finite; 
the corresponding tree amplitudes vanish making representation Eq. flFTl) senseless. For 
those finite amplitudes explicit results are presented below. Note however that in this case 
results depend on phase conventions adopted for gluon polarization vectors. To get rid 
of this dependence, for finite amplitudes we compare absolute values \An,i\ to the results 
available in the literature. All results for scattering amplitudes reported below are given in 
the FDH scheme; results in the HV scheme can be obtained using Eq. flSUI) . 

Finally, we note that QCD sub-amplitudes are often calculated using supersymmetric 
decomposition since in this way analytic computations are simplified. In particular, the 
only part of the gluon scattering amplitude that cannot be obtained from four-dimensional 
unitarity is the part where gluons scatter through a loop with virtual scalars. Dealing with 
virtual scalars is easier than with virtual gluons since the number of degrees of freedom is 
smaller. However, since our goal is to demonstrate the vitality of the method, we do not 
employ this simplification and compute sub-amplitudes, without using the supersymmetric 
decomposition. 
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A. Four gluon scattering 

We consider color-ordered sub-amplitude v4.4^i(l''*^, . . . , 4^") for the following choice of 
external momenta (p = {E,Px,Py,Pz)) 



Pi = (l,0,0,l), P2 = (1,0,0,-1) 



P3 = (— 1, sin 6', 0, cos^) , p4 



■1, — sin 6, 0, — cos( 



with 6 = 7t/3. 



Ai, A2, As, A4 


/\cut 


/\rat 


A 


- - + + 


2.53627 


0.2222 


2.75849 


- + - + 


1.90292-3.29626 i 


0.66667 


2.56959-3.29626 i 



(52) 



TABLE I: Finite parts of singular four-gluon scattering amplitudes for various gluon helicities. 
Cut-constructible and rational parts are shown separately. 



There are four sub-amplitudes that we have to consider -f- + ++, — h -|--|-, 1-+ and 

— I h. The amplitudes + + ++ and — h ++ are finite, i.e. the cut-constructible parts of 

these amplitudes vanish identically and the entire results are due to rational parts. We find 



A,i(l+,2+,3+,4+) = 2 cr X 0.33333, 
A,i(l",2+,3+,4+) = z cr X 0.7500. 



(53) 



Finite parts for the two divergent amplitudes 1-+ and — I h are given in Table [B 

The results for scattering amplitudes in Eqs. (1511531) and Table [T] are in complete agreement 



with Ref. 



4l|. 



B. Five-gluon scattering 



We consider sub-amplitude ^5 1(1^^^, . . . , 5^^^) and choose external momenta to be 

= (1,0,0,1), p2 = (1,0, 0,-1), 

P3 = e(-i,i,o,o), P4 = e(-v^, 0,1,1 

P5 = -Pl -P2-P3- P4, 



(54) 
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with ^ = 2/(1 + ^2 + 73) = 0.4823619098. We have to consider four sub-amphtudes 

+ + + + +, — 1- + + +, 1- + + and — I h +. The first two amphtudes are finite and, 

hence, are entirely due to rational parts. 



Ai, A2, . . . , As 






A 


-- + + + 


7.3382230-2.1860355 i 


0.24488559-1.4089423 i 


7.58310859-3.5949778 i 


- + - + + 


12.059206+1.7853279 i 


-7.5603579-8.4763597 i 


4.4988481-6.6910318 i 



TABLE II: Finite parts of singular five-gluon scattering amplitudes for various gluon helicities. 
The cut-constructible and rational parts are shown separately. 



For the two finite amplitudes we obtain 



2+, 3+, 4+, 5+) = i cr X (0.01056 - 0.6614 i) , 
A,i(l", 2+, 3+, 4+, 5+) = i cr X (-0.6773 - 0.4976 i) . 



(55) 



The results for the two divergent amplitudes are given in Table HIl For all scattering ampli- 
tudes, we find complete agreement with the results reported in Ref. 40|. 



C. Six gluon scattering 

We now present the results for the six-gluon scattering. We consider color-ordered sub- 
amplitude .46,i(l, 2, 3, 4, 5) for the same external momenta as considered in Ref. [l^^ Specif- 
ically, we have 

Pi = 3 (—1, cos 6' sin 0, sin6', cos 9 cos0) , P2 = —3 (1, cos 6' sin 0, sin 6', cos 6', cos0) ; 
P3 = 2 (1, 0, 1,0), p4 = ^ (1, sin/3, cos/3, 0) , 

Ps = (1, cos a sin /3, cos a cos j3, sin a) , Pe = —pi — P2 — Pz — Pa — Ps? (56) 

where 6* = 7r/4, </> = 7r/6, a = 7r/3, cos [5 = —7/19. 

We compute sub-amplitude ^6,i(l' 2, 3, 4, 5, 6) for the following helicity configurations 

+ + + + ++, - + + + ++, -- + + ++, - + - + ++,- + + -++, + ++, -- + -++, 

^ Our conventions differ from Ref. in that x and y axes are interchanged. 
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Ai, A2, . . . , Ae 


/\cut 


/\rat 


A 


-- + + + + 


-19.481065+78.147162 i 


28.508591-74.507275 i 


9.027526+3.639887 i 


- + - + + + 


-241.10930+27.176200 i 


250.27357-25.695269 i 


9.164272+1.480930 i 


- + + - + + 


5.4801516-12.433657 i 


0.19703574+0.25452928 i 


5.677187-12.179127 i 


- - - + + + 


15.478408-2.7380153 i 


2.2486654+1.0766607 i 


17.727073-1.661354 i 


-- + - + + 


-339.15056-328.58047 i 


348.65907+336.44983 i 


9.508509+7.869351 i 


-+-+-+ 


31.947346+507.44665 i 


-17.430910-510.42171 i 


14.516436-2.975062 i 



TABLE III: Finite parts of singular six-gluon scattering amplitudes for various gluon helicities. 
The cut-constructible and rational parts are shown separately. 



— I 1 1 h. Among eight amplitudes, two amplitudes + + + + ++ and — h + + ++ 

are finite; the corresponding tree amplitudes vanish. For these amplitudes we find 

2+, 3+, 4+, 5+, 6+) = i cr X (0.4968861 - 0.1838453 i) , 
A,i(l",2+,3+,4+,5+,6+) =icr X (1.0182879 + 3.0968492 i) . (57) 

The results for the finite parts of remaining six gluon scattering amplitudes are given in 
Table llllt where cut-constructible and rational parts are shown separately. For all six-gluon 
amplitudes, we find an agreement^, at least through seven digits, with the results reported 
in Ref. M. 



VII. CONCLUSIONS 



In this paper we describe a novel method for calculating one-loop scattering amplitudes, 
including their rational parts. It is based on unitarity cuts in higher-dimensional space- 
time. Similar to four-dimensional unitarity, one-loop amplitudes are obtained from tree 
amplitudes. These tree amplitudes can be efficiently calculated using recursive algorithms 
of polynomial complexity leading to an efficient method for semi-numerical evaluation of one- 
loop amplitudes. Because the method is built around integer-dimensional unitarity cuts, we 
do not foresee any difficulty with its application to chiral gauge theories. 



2 



To obtain numerical results reported Ref. [14| from our results, one should expand fi"'^ included in the 



normalization factor cr, Eq. in powers of e and substitute fi — 6. 
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Our method solves the outstanding problem of an efficient semi-numerical evaluation of 
the so-called rational part of one-loop amplitudes, an important step towards automated 
computation of NLO cross sections. The generality of the method allows a straightforward 
calculation of NLO corrections to multi-particle processes that involve virtual particles of 
arbitrary spins and masses. We hope that further development of this method will finally 
bring within reach NLO computations for such complicated processes as PP — > tt + 2 jets 
and PP ^ y + 3, 4 jets. 
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